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1.  Introduction 


One  of  the  most  important  applications  for  amorphous  polymers  is  in  transparent  armor 
packages.  Typically,  these  systems  consist  of  alternating  glass  and  transparent  polymer  layers. 
The  glass  serves  to  deform  and  fracture  the  threat  while  the  polymer  layer  acts  as  a  spall  shield 
( 1 ,  2).  Polycarbonate  (PC)  is  a  material  that  has  high  ductility  and  low  density,  making  it  a  good 
candidate  for  use  in  these  packages.  PC  is  also  being  considered  for  several  specialized 
applications  because  it  has  the  potential  to  offer  significant  weight  savings  and  better  protection 
against  relevant  threats.  To  support  ongoing  research  in  these  areas,  experiments  were 
conducted  on  PC  to  characterize  its  mechanical  response. 

The  experiments  performed  as  part  of  this  study  were  used  to  detennine  material  constants  for 
constitutive  models.  These  include  uniaxial  stress  compression  at  both  low  and  high  rates  as 
well  as  elevated  temperatures.  Low  rate  experiments  (<l/s)  were  performed  with  a  servo- 
hydraulic  load  frame  and  high  rate  tests  (~10  /s)  with  the  Split  Hopkinson  Pressure  Bar  (SHPB) 
or  Kolsky  Bar  technique.  These  experiments  were  supplemented  with  low  rate  tensile  tests. 
Material  constants  were  calculated  for  the  Johnson  Cook  (JC)  strength  model  and  the  Zerilli 
Annstrong  (ZA)  polymer  constitutive  model  (3-5). 


2.  Material 


The  specimens  used  in  this  study  were  machined  from  a  12.7-mm  thick  Bayer  MAKROLON 
General  Purpose  cast  polycarbonate  plate.  Compression  experiments  performed  on  specimens 
taken  from  the  thickness  direction  and  two  orthogonal  directions  revealed  no  evidence  of 
anisotropy. 

Table  1  lists  the  mechanical  properties  of  the  material.  The  densities  of  six  specimens  were 
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measured  using  a  water  buoyancy  method  and  the  average  value  was  found  to  be  1 197.75  kg/m 
with  a  standard  deviation  of  0. 10  kg/m  .  Longitudinal  and  shear  wave  speeds  have  been  reported 
by  Millet  and  Bourne  (6);  assuming  linear  elastic  behavior,  elastic  constants  can  be  calculated. 
The  wave  speeds  in  the  table  are  in  agreement  with  the  work  of  Carter  and  Marsh  (7)  who  report 
longitudinal  and  shear  wave  speeds  of  2. 19  km/s  and  0.89  km/s  respectively. 
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Table  1.  Material  properties  for  PC.  aDenotes  value  calculated 
by  Millet  and  Bourne  (6). 


Density  (kg/m3)= 

1197.75 

C  Longitudinal  (kl7l/s)  — 

2.13 

aCShear  (km/s)= 

0.88 

Modulus  of  Elasticity  (GPa)= 

2.59 

Shear  ModulusfGPaj^ 

0.93 

Lame  Constant  (GPa)= 

3.58 

Bulk  Modulus  (GPa)= 

4.20 

3.  Experimental  Procedures 


3.1  Low  Rate  Compression  and  Tension 

Quasi  static  compression  tests  were  conducted  using  an  Instron  model  no.  1331  servo-hydraulic 
load  frame.  During  the  tests  a  22.2  kN  load  cell  measured  the  force  applied  to  the  specimen  and 
a  linear  variable  differential  transformer  (LVDT)  measured  the  crosshead  displacement.  *  Strain 
rates  using  this  setup  varied  from  0.005/s  up  to  0.4/s.  The  specimens  were  cylindrical  and  had 
nominal  dimensions  of  5X5  mm  (lengthxdiameter).  The  PC  samples  were  loaded  between  two 
tungsten  carbide  platens  and  the  contact  surfaces  were  lubricated  with  M0S2  grease  to  reduce 
friction.  To  perfonn  low  rate  tests  at  elevated  temperatures  an  Instron  SFL  Temperature  Control 
Chamber  was  used.  Specimens  were  compressed  at  temperatures  of  38  °C,  54  °C,  and  71  °C. 
During  temperature  tests  the  platens  and  actuator  anns  were  heated  to  the  desired  temperature 
before  the  specimen  was  positioned.  Once  in  place,  the  sample  was  allowed  to  rest  for  15  min 
before  testing  to  allow  the  temperature  to  equilibrate. 

Low  rate  tension  experiments  were  completed  using  the  same  Instron  model  no.  1331  machine 
previously  described.  A  dimensioned  drawing  of  the  tension  specimen  can  be  found  in  figure  1 . 
To  accommodate  these  specimens  threaded  collars  were  attached  to  actuator  anns  of  the  load 
frame. 


Machine  compliance  corrections  were  incorporated  into  all  data  analysis. 
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A  Retiga  2000-R  digital  camera  was  used  to  record  the  deformation  of  the  samples  during  tests 
and  to  provide  a  confirmation  of  measured  strain. 

3.2  High  Rate  Compression 

High  rate  compression  experiments  were  conducted  using  a  SHPB  (8-10).  For  strain  rate 
experiments  up  to  ~6500/s,  7075-T6  aluminum  bars  with  diameters  of  9.525  mm  were  used.  In 
order  achieve  even  higher  rates  (-15, 000/s)  specimens  with  a  reduced  gage  length  were  tested 
using  bars  with  a  diameter  of  4.76  mm.  Cylindrical  samples  were  tested  in  both  setups  and  the 
associated  dimensions  were  5X5  mm  and  2x2  mm  respectively.  These  dimensions  were  chosen 
so  that  an  aspect  ratio  of  1  would  be  maintained  for  all  compression  tests.  Each  bar  was 
instrumented  with  two,  diametrically  opposed  strain  gages'!  to  measure  the  elastic  waves  within 
the  bars.  As  with  the  low  rate  tests,  the  contact  surfaces  of  the  PC  samples  were  lubricated  with 
M0S2  grease.  The  data  analysis  for  these  experiments  includes  a  correction  for  wave  dispersion 
in  the  bars  (11—14).  During  all  dynamic  tests  adiabatic  behavior  is  assumed. 

Dynamic  compression  tests  were  also  conducted  at  elevated  initial  temperatures  to  investigate 
the  combined  effects  of  strain  rate  hardening  and  thermal  softening.  An  insulated  chamber  was 
constructed  that  enclosed  the  specimen  and  -100  mm  of  both  the  incident  and  transmission  bars. 
Hot  air  was  then  circulated  throughout  the  chamber  for  -20  min.  Because  soldering  a 
thermocouple  directly  to  the  PC  specimen  was  not  possible,  one  was  placed  within  10  mm  of  the 
sample  to  measure  ambient  air  temperature.  It  is  assumed  that  given  the  soak  time  the  air 
temperature  at  this  location  provides  a  reasonable  approximation  of  specimen  temperature. 


1  Two  series  of  Vishay  Micromeasurement  strain  gages  were  used:  WK-06-250BF-10C  and  EA-13-062AQ-350. 
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4.  Results 


4.1  Low  Rate  and  High  Rate  Compression  at  Room  Temperature 

Quasi  static  experiments  were  conducted  at  three  strain  rates:  0.005/s,  0.03/s,  and  0.4/s.  Figure  2 
shows  a  representative  stress  strain  curve  at  each  rate.  All  three  stress-strain  curves  in  the  figure 
display  similar  behaviors.  First,  stress  increases  nonlinearly  up  to  a  peak  strain  of 
-0.08  mm/mm.  It  is  assumed  that  this  behavior  is  visco-elastic.  Beyond  this  peak  the  flow  stress 
begins  to  decrease  as  a  result  of  strain  softening.  At  a  strain  of  0.35  mm/mm,  the  flow  stress 
begins  to  increase  again  due  to  strain  hardening.  These  behaviors  occur  because  of  molecular 
chain  rotations  and  translations  and  have  been  discussed  in  detail  by  Mehta  and  Prakash  (75). 
This  data  suggests  that  the  yield  stress  and  to  a  lesser  extent  the  elastic  modulus  (as  indicated  by 
the  assumed  elastic  loading)  are  positively  correlated  to  strain  rate.  Specimens  were  inspected  at 
the  conclusion  of  each  test  and  no  evidence  of  failure  could  be  found.  The  results  presented  here 
are  in  agreement  with  the  work  of  previous  authors  (16-18). 


True  Stress  v.  True  Strain 


True  Strain  (mm/mm) 


Figure  2.  Representative  stress  -  strain  curves  from  low  rate  experiments. 

SHPB  compression  tests  were  done  at  four  different  rates  and  figure  3  shows  representative 
results  for  the  tests.  The  same  three  behaviors  noted  in  the  low  rate  experiments  are  found  here. 
The  test  conducted  at  1750/s  does  not  show  evidence  of  strain  hardening,  but  this  is  because  the 
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specimen  was  only  loaded  to  a  strain  of  ~0.3  mm/mm.  Also,  the  stress-strain  curves  did  not 
return  to  zero  stress  because  only  a  portion  of  the  specimen  unloading  was  recorded.  Recovered 
specimens  were  visually  inspected  and  there  was  no  evidence  of  fracture.  At  some  point 
between  quasi-static  and  dynamic  conditions,  the  strain  rate  effect  on  the  apparent  yield  stress 
becomes  decidedly  more  dramatic.  This  phenomenon  occurs  due  to  the  activation  of  a  secondary 
molecular  process  at  high  strain  rates  (19-20).  Consequently,  the  polymer  chains  become  stiffer 
and  higher  stresses  are  achieved. 


True  Stress  v.  True  Strain 


True  Strain  (mm/mm) 


Figure  3.  Representative  stress  -  strain  curves  from  SHPB  compression  tests. 

In  figure  4  the  observed  yield  stress  of  PC  is  plotted  as  a  function  of  the  strain  rate  to  determine 
an  approximate  activation  strain  rate  for  the  secondary  process.  Logarithmic  trendlines  were 
fitted  to  the  experimental  data  and  extrapolated.  The  intersection  of  the  extrapolated  lines 
provides  a  rough  estimate  of  where  the  transition  to  the  more  pronounced  rate  effect  occurs. 
Based  on  the  data  the  threshold  strain  rate  is  — 10/s,  but  others  have  reported  it  closer  to  100/s 
(20).  Additionally,  it  appears  that  beyond  3500/s,  higher  strain  rates  begin  to  have  a  diminishing 
effect  on  the  yield  behavior. 
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Figure  4.  Yield  stress  as  a  function  of  strain  rate.  Extrapolation  of  linear  fits  (dotted)  suggest  a 
transition  to  enhanced  hardening  at  ~10/s. 

4.2  Compression  at  Elevated  Temperatures 

Select  results  from  elevated  temperature  tests  are  shown  in  figure  5.  Predictably,  higher 
temperatures  result  in  decreased  yield  and  flow  stresses.  Figure  6  shows  the  yield  stress  from  all 
the  experiments  as  a  function  of  temperature.  Low  temperature  tests  (i.e.,  below  room 
temperature)  were  not  conducted  as  part  of  this  program,  but  several  researchers  have  performed 
such  experiments  (16,  18,  19).  They  found  that  lower  temperatures  have  the  reverse  effect,  i.e., 
increased  yield  and  flow  stresses.  Similar  to  the  high  rate  tests,  these  increases  are  attributed  to 
the  restriction  of  secondary  molecular  motions. 
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Figure  5.  Results  from  compression  tests  at  elevated  temperatures. 


Figure  6.  Yield  stress  as  a  function  of  temperature. 
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4.3  Low  and  High  Rate  Tension 

Figure  7  shows  engineering  stress-strain  curves  from  low  rate  tensile  tests  using  the  dog  bone 
specimen  shown  in  figure  1 .  A  singular  compression  test  is  included  in  the  figure  for 
comparison.  As  in  compression,  the  tensile  data  shows  regions  of  visco-elastic  loading, 
softening,  and  hardening.  An  interesting  observation  from  the  experiments  was  the  tendency  of 
the  neck  to  propagate  along  the  gage  length  as  shown  in  figure  8.  The  movement  of  the  neck 
along  the  length  of  the  specimen  corresponds  to  the  relatively  flat  portion  of  the  stress  strain 
curve  between  0.1  mm/mm  and  0.4  mm/mm.  Once  propagation  has  finished  the  flow  stress 
begins  to  increase  due  to  hardening.  A  detailed  discussion  of  this  behavior  has  been  presented 
by  Wu  and  Van  Der  Giessen  (21). 


Eng.  Stress  v.  Eng  Strain 


Figure  7.  Results  from  low  rate  tension  experiments. 

Dynamic  tension  experiments  were  not  performed  as  part  of  this  study,  but  the  high  rate  tensile 
response  of  PC  has  been  documented.  Cao  et  al.  (22)  conducted  high  rate  tension  experiments 
and  found  that  tensile  yield  and  flow  stresses  are  positively  correlated  to  strain  rate.  They  report 
that  for  tests  at  rates  between  370/s  and  1700/s  yield  tends  to  occur  at  roughly  10%  engineering 
strain  at  a  stress  of  100  MPa.  Sarva  and  Boyce  (23)  also  performed  high  rate  tension 
experiments  on  PC.  They  were  able  to  initiate  two  necks  in  straight  gage  specimens  by  altering 
the  diameter  of  specimens  while  keeping  the  gage  length  constant.  Additionally,  they  report  that 
the  formation  of  necks  depends  on  the  loading  velocity  and  method  used  to  grip  the  specimen. 
Finally,  Sarva  and  Boyce  successfully  simulated  the  results  from  their  experiments  using  the 
ABAQUS  finite  element  code. 
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Figure  8.  Neck  propagation  in  straight  tension  specimens. 


5.  Constitutive  Modeling 


The  experimental  data  was  fit  to  the  JC  and  ZA  (polymer)  constitutive  equations.  The  JC 
equation  is  written  as: 

a  =  (A  + Bspn)(l  + (1) 


where  T*  is  the  homologous  temperature: 

rp*  _  T~TR 

Tm~T  r 


(2) 


In  the  equation  A,  B,  C,  m,  and  n  are  all  constants  related  to  the  material.  T  is  the  material 
temperature  and  TR  is  a  reference  temperature,  which  in  the  present  case  is  equivalent  to  room 
temperature.  Tm  is  usually  the  melting  temperature;  here  it  is  used  as  free-parameter  to  improve 
the  quality  of  the  fit.  e  is  the  strain  rate  and  £0  is  the  reference  strain  rate  for  which  1/s  was  used. 
This  form  of  the  JC  constitutive  model  is  readily  available  in  most  hydrocodes  and  for  that 
reason  it  is  often  chosen  to  model  material  behavior. 


To  determine  the  material  constants  a  series  of  least  square  fits  to  the  experimental  data  are 
required.  During  this  process,  only  the  compression  data  is  used.  Assuming  T  =  Tr  =  295  K  and 
sp  ~  0  values  for  A  and  C  can  be  estimated  from  a  logarithmic  fit  to  the  data  in  figure  4. 

Estimates  of  m  and  Tm  can  be  obtained  from  a  linear  fit  to  the  temperature  data  in  figure  6.  In 
this  case  a  fit  to  the  low  rate  data  was  used  because  a  wider  range  of  temperatures  were  available. 
As  such  it  was  assumed  that  e=0. 005/s  and  again  that  sp~0.  The  remaining  terms  B  and  n  can  be 
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fit  to  the  stress-strain  curves  in  figure  3  assuming  isothennal  behavior  at  low  rates  (0.005/s- 
0.4/s)  and  adiabatic  conditions  at  high  rates  (>0.4/s).  For  the  adiabatic  curves  the  specific  heat 
capacity,  Cv,  was  taken  to  be  1 .3  KJ/kg-K.  The  ratio  of  plastic  work  converted  to  heat,  P,  was 
assumed  to  be  0.5  based  on  the  work  of  Li  and  Lambros  (24).  Because  the  majority  of  the 
constants  were  estimated  from  yield  data  it  is  necessary  to  adjust  each  one  to  that  value  that  best 
represents  the  experimental  data. 

Table  2  lists  the  JC  parameters  and  figure  9  shows  JC  stress  strain-curves  compared  to  the 
experimental  data.  The  JC  constitutive  equation  was  not  intended  to  capture  the  softening  that 
occurs  after  specimen  yield  and  as  such  it  does  not  represent  the  low  rate  (isothermal)  data  well. 
However,  at  high  higher  rates  where  adiabatic  conditions  are  assumed  the  equation  provides  a 
reasonable  approximation  of  the  experimental  data.  This  is  because  the  introduction  of  thennal 
softening  into  the  equation  diminishes  the  strain  hardening  effects  resulting  in  a  better 
approximation  of  the  constitutive  behavior.  Note  that  while  activation  of  thennal  softening  in 
this  equation  provides  a  better  approximation  of  the  data,  this  is  not  the  correct  real  world 
mechanism. 


Table  2.  JC  parameters. 


A= 

80  MPa 

B= 

75  MPa 

C= 

0.052001  — 

m= 

0.548  — 

n= 

2  — 

Tmcit= 

562  K 

P= 

0.5  — 

P= 

1220  kg/m3 

cv= 

1.3  KJ/(kg  K) 

The  experimental  data  was  also  fit  to  a  simplified  form  of  the  ZA  polymer  constitutive  model 
that  assumes  no  pressure  dependence.  While  it  is  known  that  PC  has  some  pressure  dependence 
there  is  not  enough  available  data  to  fit  that  aspect  of  the  model.  To  estimate  parameters  for  the 
simplified  model  the  method  outlined  by  Casern  (25)  was  followed.  Table  3  lists  the  constants 
for  the  ZA  polymer  constitutive  model  and  figure  10  compares  the  ZA  stress-strain  curves  to  the 
experimental  data.  The  model  reasonably  captures  the  material  behavior,  but  significant 
variation  persists  despite  the  additional  parameters.  It  was  possible  to  adjust  the  constants  so  that 
they  provided  a  very  accurate  representation  for  either  the  assumed  isothermal  or  adiabatic 
conditions,  but  not  both.  The  inability  of  both  models  to  represent  material  behavior  over  the 
entire  range  of  strain  rates  occurs  because  there  is  no  mechanism  in  the  equations  to  account  for 
the  transition  to  enhanced  rate  hardening  as  a  result  of  the  secondary  molecular  process. 


10 


True  Stress  v.  True  Strain 


True  Stress  (MPa) 


Figure  9.  Comparison  of  JC  calculations  to  experimental  data. 


Table  3.  ZA  model  parameters. 


B0= 

0.006715948 

1/K 

Bi= 

0.00009503 

1/K 

BDa= 

550 

MPa 

-^opa- 

48 

MPa 

®a= 

-8 

— 

CUb= 

-0.01 

— 

P= 

0.5 

— 

a0= 

0.00655 

1/K 

a1= 

0.00004 

1/K 

P= 

1220 

kg/rn 

cv= 

1.3 

KJ/(kg  K) 
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True  Stress  v.  True  Strain 


True  Strain  (mm/mm) 


Figure  10.  Comparison  of  ZA  calculations  to  experimental  data. 


6.  Discussion  and  Conclusion 


A  series  of  mechanical  tests  were  conducted  on  polycarbonate  to  characterize  its  mechanical 
response.  These  included  experiments  in  tension  and  compression  and  covered  a  strain-rate 
range  of  0.001/s  to  15k/s  and  a  temperature  range  from  295  to  343  K.  It  was  shown  that  strain 
rate  and  temperature  are  important  factors  that  influence  the  material  behavior.  Data  from  the 
experiments  were  used  to  calculate  parameters  for  the  JC  strength  model  and  a  simplified  ZA 
polymer  strength  model.  Both  models  were  able  to  capture  the  basic  features  of  the  mechanical 
response  and  should  be  suitable  for  certain  classes  of  impact  problems. 

Some  significant  drawbacks  of  the  work  exist.  First  and  foremost  is  the  absence  of  high-pressure 
data.  The  yield  behavior  of  polycarbonate  has  been  shown  to  be  pressure  dependent.  This  has 
been  established  from  plate  impact  experiments,  e.g.,  Millet  and  Bourne  (6),  who  found  that  the 
shear  strength  of  the  material  increases  with  pressure  under  uniaxial  strain  shock  loadings. 
Limited  work  has  been  performed  with  the  pressure  shear  plate  impact  method  (PSPI)  (26-28). 
Mehta  and  Prakash  (15)  have  performed  PSPI  experiments  on  PC  and  found  that  a  pressure  of 
1  GPa  and  strain  rate  of  400,000/s  result  in  a  significant  increase  in  the  yield  (0.1  GPa)  and  flow 
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stress  of  the  material.  Furthermore,  the  authors  reported  that  the  material  does  not  undergo  any 
strain  softening  under  these  conditions. 

In  addition  to  the  added  component  of  pressure,  these  experiments  have  the  further  benefit  of 
increased  strain-rates  (beyond  15  k/s).  The  higher  rate  data  are  beneficial  for  certain  classes  of 
ballistic  problems,  and  is  another  shortcoming  of  the  current  work.  Other  methods  to  achieve 
high  rate  that  would  be  appropriate  for  this  material  are  miniature  Kolsky  bar  methods,  for 
example  Casern  et  al.  (29)  have  designed  a  miniaturized  SHPB  that  allows  experiments  to  be 
conducted  at  strain  rates  in  excess  of  150,000/s.  Triaxial  stress  experiments  with  both  servo- 
hydraulic  and  SHPB  methods  can  also  be  employed  to  study  the  effects  of  pressure  over  a  range 
of  strain-rates. 

In  future  work  it  is  suggested  that  a  more  thorough  constitutive  model,  such  as  the  one  presented 
by  Mulliken  and  Boyce  (30),  be  used  to  represent  the  data.  Finally,  experiments  should  be 
conducted  to  describe  the  fracture  behavior  of  polycarbonate  in  a  ballistic  environment. 
Although  more  appropriate  for  metals,  the  JC  fracture  model  could  offer  a  place  to  start. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


DSI 

Dynamic  Science,  Inc. 

JC 

Johnson  Cook 

LVDT 

linear  variable  differential  transformer 

PC 

Polycarbonate 

PSPI 

pressure  shear  plate  impact 

SHPB 

Split  Hopkinson  Pressure  Bar 

ZA 

Zerilli  Armstrong 
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